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ABSTRACT 


An analysis of Synchronous Time Division Multiplexing is presented. Packets 
of information arrive at the system as a compound Poisson process, and can 
be transmitted only during individual periodic intervals. Packet arrivals 
may be blocked (lost) if the system has a finite capacity and is congested. 
Using the theory of semi -regenerative processes, the distribution of the 
number of packets in the system (system size) is found. This nonstationary 
distribution is used to determine the complete system behavior, including 
the delay distributions, the blocking probability, and the density of the 
system size at arrival instants. Numerical examples illustrate applications 
of the results given. 


L I TO DUCnO H 


Spacecraft data Is collected by means of a variety of sensors which typically 
operate simultaneously. The data from the various sensors must be funneled 
through a common telemetry channel . The current method of combining the 
data is by Synchronous Time Division Multiplexing (STDM). 

This paper represents the initial phase of research in the study of Deep 
Space packet telemetry techniques [1], This study, which is part of the 
NEEDS (NASA End-to-End Data System) program, will identify and compare 
various strategies for managing, packetizing and multiplexing spacecraft 
data. The baseline technique defined in [1] is STDM. 

A unified approach to analyze the performance of STDM is given herein. 

Using this approach, several new results are obtained for a large cliiss of 
arrival streams, particularly for the finite buffer capacity case. 

The design engineer may use these results to completely predict system 
behavior, including channel utilization, probability of data loss due to 
buffer overflow, and queueing delay. By varying the design parameters of 
the model, the appropriate channel allocation and buffer sizing can be 
determined such that the performance will meet the prescribed criteria for 
each source. In addition, sensitivity analyses can be done by exploiting 
the generality of the arrival process model. 

In STDM, each source is assigned a fixed sequence of time intervals during 
which it may transmit information. The assignment is predetermined and 
does not adapt to fluctuations in the traffic load. In the definition of 
STDM in this study, a structure is imposed on the time intervals assigned 
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to the sources, and fs defined as follows. Time is divided into equal 
length frames. Each source is allotted a fixed contiguous portion of 
each frame, defined as a slot, which remains in a fixed ordered position 
from frame to frame. While the slot sizes are constant for a given source, 
they may vary between sources. The amount of Information transmitted in 
a slot is called a packet (for that source). 

It is clear that STDM provides excellent performance if the data arrival 
times are deterministic, one packet per frame, In the past, spacecraft 
data were collected in a predetermined fashion so that STDM was well suited 
for the environment. However, future missions will use source encoders 
as well as sensors that are triggered by random-time events. These will 
cause the information packets to arrive at random times, in which case 
STDM may not be as efficient as other alternatives. This paper investi- 
gates the performance of STDM with random time-of-arrival input streams. 

In any multiplexing scheme the source data are colocated at the concentrator. 
This allows the possibility of sharing another finite resource; buffer 
space. However, in this paper we will assume each source has its own 
individual buffer. This important assumption implies that the queueing 
behavior of any source is independent of all other sources, and only one 
source needs to be considered. Since the STDM scheme considered here does 
not utilize any of the advantages due to the colocation of the sources, it 
is also referred to as Time Division Multiple Access (TDMA). 

The analysis of the STDM queueing process presented here is based on the 
theory of semi -regenerative processes, assuming a Poisson message arrival 
stream. Messages will consist of groups of packets, where the distribution 


3 


of the group size Is arbitrary (i.e. a compound Poisson packet arrival stream 
Is assumed). The quantities of Interest are: the number of packets or 
messages In the system ("system size"), the total time spent In the system 
("delay"), and In the case of a finite capacity system, the probability of 
being blocked. 

In this paper we give new results describing STDM/TDMA behavior, Including 
the actual delay distributions and a complete characterization of Intraframe 
system behavior for finite packet capacity systems. A list of previous work 
on STDM/TDMA Is given by [2]-[8] and [13]-[17]. Birdsall et al. [13] and Dor 
[2] found the system size distribution at framing Instants as well as 
the blocking probability for a finite capacity system. Chu [14] extended 
these results to the case of compound Poisson arrivals. Konhelm [3]. Hayes 
[4]. Lam [5] and Rubin [6.7] all concentrated on the infinite buffer capacity 
case. Konheim [3] derived the packet system size distribution at framing 
instants and the mean virtual packet waiting time. Hayes [4] describes intra- 
frame behavior, and gives the actual waiting time distribution as well as the 
packet system size distribution. The first work to consider the packet trans- 
mission time within the frame was Lam's [5]; his paper gives the message system 
size distribution and mean message delay. Van [8] extended Lam's work to allow 
a finite message capacity. Rubin [6,7] used a discrete-time model (with infi- 
nite buffer capacity) to obtain the system size and delay distributions. His 
model is more general than the one used here in that a source is allowed multi- 
ple contiguous slots in each frame. In addition, Rubin points out that the 
interarrival times may be correlated within the discrete time unit, so that a 
more general arrival process is allowed. Generalizations in other directions 
are given by [15]-[17]. 
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The STDM model will be defined In the next section. Section III will 
present the fundamental analysis leading to the key result (28)-(29), the 
steady state packet system size density. This is used to find the packet 
system size at arrival instants (a non-Markov chain), as well as the 
blocking and truncation probabilities. Section IV contains the deriva- 
tion of the packet and message delay distributions. Transform relations 
are given in Section V for the infinite capacity case. Numerical examples 
are then given in Section VI, followed by a concluding summary. The 
appendix contains an analysis of the message system size for a possibly 
finite message capacity system; this extends the works of Van [8] and Lam [5]. 
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II. DEFINITION OF THE STDM MODEL 


In STW1, time Is divided Into equal length frames » Each frame Is further 
divided Into ordered slots, and one ? lot Is assigned to eacb source. Let 

T * frame duration 

M = number of sources being multiplexed 
=» proportion of frame assigned to source m, m » 1, 2, M 

so that 

« slot duration for the mth source, m = 1, 2, M 

The STDM structure is illustrated In Figure 1. Frequently, STDM (or TDMA) 
Is used with all the sourcfe^ viavlng the same slot size, so that 6^^ = 1/M, 

As was previously mentioned, the queueing behavior for a given source is 
independent of all the other sources, so that only one needs to be con- 
sidered. Define the time origin as though a slot for the source under 
consideration has just ended. Thus, the source may use the cross-hatched 
slots of duration 6j^T illustrated in Figure 1. 

All subsequent definitions will refer only to the individual source under 
consideration. For example, the "packet system size" is the total number 
of packets in the system belonging to that source. In addition, the 
subscript on Sj^Will be dropped, so that 

6 = proportion of frame assigned to the source 


Messages enter the system as a Poisson arrival process. Thus, letting 

Ajj * time of th{* nth message arrival 
X » message arrival rate, messages/time unit 

we have that the Interarrival times are Independent and have the distribution 

P(\+l “ j< t) * 1 - e"^* , t ^ 0, any n 

Messages are assumed to consist of a positive Integer number of packets. 

Let 


® number of packets In the nth message 

Then {6^,} Is a sequence of Independent and Identically distributed random 
variables with finite mean and an otherwise arbitrary density given by 

“ P ~ , A®1,2, ...,n’‘l,2, ... 

A packet is transmitted whenever the system is not empty at the beginning 
of a slot. If a packet arrives to an empty system during a slot interval, 
it must wait for the next slot before it can begin being transmitted (I.e. 
a complete slot is required for transmission)* 

The "system size” can be counted either in packets or in messages, and is 
defined to include both those in the buffer and that "in service" (i.e. 
being transmitted). Figure 2a illustrates a packet system size sample 
path, while Figure 2b shows the message system size for the same arrival 
stream. The sequences {R^} and represent the departure instants of 
packets and messages respectively. Each of the systems illustrated has a 
finite capacity. Packets or messages that are blocked are assumed to be 
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Irretrievably lost. From Figure 2 it is clear that placing a finite 
constraint on the packet state space implies a different blocking procedure 
for the system than placing a finite constraint on the message state space. 

The primary concern of this paper is the determination of the packet system 
size distribution (although the message system size distribution is derived 
in the appendix). Therefore we define for the body of this paper 

X|. * packet system size at time t 

N = capacity of the system in packets (possibly infinite) 

For N < « the implicit blocking procedure is as follows; if a message of 
G packets arrives to a system containing N-L packets, and L is less than 
6, then L packets will be accepted and the remaining G-L will be blocked 
(lost). Thus message integrity is not maintained, and several partial 
messages could be in the buffer at the same time. (The term "message" 
seems to imply an underlying mutual information between packets in the 
same arrival group, which may not be the case. The packets may actually 
contain independent information, with multipacket "messages" being used 
solely to model the burstiness of the arrival process.) 

The next section begins by defining processes which are "embedded" within 
the {X^} process. These processes are then characterized, and provide the 
means for determining the distribution of X^. 


III. PACKET SYSTEM SIZE ANALYSIS 


In this section we utilize the theory of semi -regenerative processes to find 
the distribution of the packet system size X^.. We will identify the proc- 
esses and functions necessary for the development of tho theory, and outline 
the approach leading to the main result ((28), (29)). Some immediate 
applications are then provided, including computation of the blocking 
probability. 


The following processes are embedded within the process {X^}: Let 


= time of the nth packet departure 

Y„ = Xd+ = packet system size just after the nth departure 
n 


and define 


Zf = Y„ for R„ , < t < R„ 

t n n-1 — “ n 

Since the arrival process is Poisson, {Y|^> is a (time-homogeneous) Markov 
Chain (MC), {(Y^^, R^^)} is a Markov Renewal Process (MRP), and {Z^} is the 
minimal Semi-Markov Process (SMP) associated with {(Y^^, R^)}. {X^} is 

a Semi -Regenerative Process (SRP) with respect to the MRP ((Y^, R^))- 
These facts can be easily checked. For "definitions of the above terms 
see [9], 


An illustrative example showing the various processes defined above is 
given in Figure 3. 
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Before proceeding further, let us first define a convenient nitatlon for 
convolutions of the density {g^} of the number of packets per message: let 


g*'^ = P(G, + Gg + • • • + Gj » a) , A * 1 , 2, . . . , J , 2, . . . 


and let 


-*o 


1 If A - 0 

0 If A I* 0 


Now consider the MC Let the transition probabilities for {Y^} be 

denoted 


" j|^n " Oli* j iN - 1 


( 1 ) 


Then for 0 < 1 £ N - 1 , 0 £ j < N - 1 


0-1+1 / k 

P|k messages arrive In [0, T] and^^ G^ = j - 1 + 1 


k=0 


n=l 


j-i+1 

Z -xT mZ. *k If j - 1 + 1 > 0 
k! *j-i+l 
k=0 


( 2 ) 


otherwise 


For 1 = 0, 0 < 0 < N - 1 we have 




j+1 

-E^r 

r=l Jo 


P(Yl = jIYq = 0. ^ G^ = r) xe'^^dt 


10 




j+1 

E«,E 

r»l m=0 

J '*T j-r+1 

E 

T-6T k=0 


'T-6T j-r+1 k 

0 k=0 




j+1 j-r+1 _*k 

r=l k=0 
- UT6)'^+^)| 

i 

The case j = N - 1 is found using 


UT)'^^’ + 


- (e-^^[xTll+6)]'^^^ 


(3) 


N-2 

^iN-1 ' ^ ■ S ^ij 

j=0 

or by summing (2) or (3) from j = N - 1 to infinity. 

Define 

P = X6T (5) 


to be the traffic intensity, where 
G = E(G^) 

is the mean number of packets in a message. By assumption G < «. 

The MC {Yj^} is irreducible and aperiodic, and is positive recurrent if and 
only if either N < « or p < 1. In this case the stationary distribution 

TT.- = lim P(Y„ = j(Yo = i) j = 0, 1. ..., N - 1 
J rH« 


for the chain exists independent of i, is unique, and satisfies 
N-1 

'H’j “ (6) 

i=0 

N-1 

L 'j ' ’ 

j=o 

If N < ®, Van [8] gives an efficient method for evaluating {ir.}. 

J 

Now consider the MRP {(Y^, R^)}. Let 

Qij(t) = = j, R^^^ - RplMYn = i) 0 1 i. j < N - 1 (8) 

be the Semi -Markov Kernel (SMK) for (Y, R). Also, for i, j fixed let 
{Rf,(ij j)} be the sequence of times at which a departure occurs and j are 
in the system immediately thereafter, given that we start with i in the 
system. This is a delayed renewal process for each i, j. Let 

N..(t) = sup {n: R»(i. j) < t} (9) 

n>0 " 

be the associated counting process to the process {R^(i, j)}. Let 

niij.(t) = E(N^.j(t)) (10) 
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be the average number of visits to state j by time t, given we start in 
state i, Then is the Markov Renewal Kernel (MRK) for the MRP (Y» R), 
The MRK can be written in terms of the SMK (useful for transient results): 

00 

mij(t) (11 

n=0 


Where 

= KUn = Vn i ‘I''. = 

is the n-stage transition distribution, which satisfies 


N-1 rt 

‘>k"j 

k=0‘^o 


(t - s)dQ^i^(s) 


( 12 ) 


For each i, j the process (R^(i, J’)) is periodic, so that each state of 
(Y, R) is periodic with period T. 

R^(j, j) is the recurrence time for state j. Note that we are speaking of 
the state space of (Y, R), which is easiest to visualize in terms of (Z^}. 
Thus Ri(j, j) is the time spends between two successive visits to 
state j. 

Let 

nj = |E[R.j(j, j)]}"’’ j = 0, 1 N-1 (13) 

be the mean recurrence rate for state j. 
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Let 


■"j ■ 


N - 1 


04) 


i » 0, 1* 


■ • » » 


be the mean sojourn time in state j, i,e. the mean time spends in 
state j before going to any other state. Then ([9], p. 329) 


* rr j ' 0, 1 1 . . . , N - 1 
m 


where 



05) 


(16) 


is the mean steady state sojourn time (for 

(It can be shown that the steady state probability density of {Z^} is 
given by {nijHj}. independent of t, even though the SMP {Z^.} is periodic, 
in our examination of the SRP {X^}. however, we will find that the 
periodicity of {X^} is evidenced by the periodicity of its steady state 
distribution.) 

Having established the necessary machinery, we may finally investigate the 
SRP {X^}. Denoting 

Pi(-) * P(*|Yo * 0 . 

we have 

P^(X^ = j) = P^(X^ = j, R, > t) + P^(X^ = j, R^ 1 t) . 
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Now 


N-1 ft 

P^(X^ » j, £ t) I P^(X^ » j|Y^ - k, » s) dP^(Y^ » k, R^ s s) 

k^Jo 


N-1 /*t 

-£ I '■k<vs = 

k =0 Jo 


j) dQ^i^Cs) 


since {X|.} Is semi regenerative. Thus {X^} satisfies a Markov Renewal 
Equation (MRE): 

N-1 rt 


P^(X^ = j) = P|(X^ = j, R^ > t) + 


N-l 

E I 

k=0 Jo 


= j) dQ^,^(s) (17) 


The solution (which can be shown is unique) to this MRE is ([9], p. 324) 


N-1 ft 

>,(Xt = j) =52 

k=0 Jo 


j, R^ > t - s) dm^,^(s) 


(18) 


We are interested in the steady state behavior of {X^}. If either N < « 
or p < 1 then the (periodic) MRP (Y, R) is irreducible and recurrent, so 
we may apply the Key Renewal Theorem ([9], p. 334) to obtain 

N-1 °» 

y; '’i(w = J) = "kE ’’ki w "'T) (19) 

“To 

where 0 _< t < T and i , j = 0, 1 , . . . , N - 1 . Using (15) and defining 

Pi^’^r+mT = • 0 < T < T, j = 0, 1, .... N (20) 

nt^ • ^ Mil J 
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we have: if N < ® or p < 1 , then for 0 < t < T 


N-1 

^ ~ XI *’k^^+mT ' > T + ml) 


( 21 ) 


k=0 m=0 


j « 0, 1, ..., N - 1 


and 


N^l 

j=0 

If k > 0, j < N then 


( 22 ) 


’’k* **1+1111 ' ■* 


while for k = 0, j < N 


ifm = 0, j-k>0 

j, > t+idT) = J 

I (23) 


otherwise 


^O^^T+mT " j> R] ^ "’T) 


_ _*A p-Xx (At)^ . p-A(t+6T) [a(t+5T)]^ \ 

■Z_i®J r *! 1 . «->T ® «! I 

A=0 ( ' ® f 


(24) 


For j > 0 the mean sojourn time is simply 




( 25 ) 
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ijmlMtr' 


itIq * T • P (at least 1 arrival In [0, T - fit)) 


00 

+ ^ (W + 1)T • P 
m«l 



(none arrive In [0, mT - 6T) and at least 
1 arrives In [mT - «T» (m + 1)T - 6T)) 


(26 


Thus 


m = T 



e-XT(l-5)) 

TT?^ 

. 


(27 


Combining (21)-(24) we obtain the primary result: If N < ® or p < 1 then 
the steady state probability density of exists, m Integer, 0 £ t < T, 
and 1s given by 



Xf^(t) = 1 



(29 


where m is given by (27). 
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The formulas (28), (29) provide the basis for the determination of the 
various quantities that describe STDM queueing behavior at steady state. 
Having found the "virtual" packet system size distribution, we next find 
the "virtual" truncation and blocking probabilities. We then complete 
the section by investigating these same quantities from the perspective 
of an "actual" arrival. 


If a message arrives at a time t = mT + t and finds X = N - J packets in 
the system (N * capacity), then the message will be truncated if it contains 
more than J packets. Our definition of truncation includes the possibility 
of J = 0, i .e., the probability of being truncated includes the probability 
X|^(t) that the entire message is lost. The steady state virtual message 
truncation density is given by 


>t(t) £ 

j=0 J!;=j+1 


g 0 < T < T 


(30) 


where the term "virtual" is used to indicate that this event Is conditioned 
on the arrival occurring at a time x + mT, m integer. 


Consider the probability that a packet is blocked, where the packet is 
chosen uniformly from all packets which attempt to enter the system. 

Suppose the packet arrives within a message of length L packets, which 
arrives at time t = x + mT and finds N-J in the system. The packet will 
be blocked if L > J and the packet is one of the L-J packets which are 
not accepted. Since the packet is randomly chosen, 

tg. 

P (packet arrives in a message of size z) = — (31) 
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and 


(32) 


P (packet is one of the (t-j) not accepted out of l) 



Theirefore, the steady state virtual packet blocking probability is 


Mt) 


N 

E 

j-o 




E 

<.sj+l 


t - .1 

5 ‘ 


0 < T < T 


(33) 


This represents an "average" virtual blocking probability, with no prior 
knowledge relating to a bias in the selection procedure. If an ordering 
exists on the packets within a message, and if the selection of which are 
blocked or not is based on the ordering, then the appropriate (nonuniform) 
distributions may be used to determine the blocking probabilities for each 
packet position. 


In Appendix B it is shown that for large n, the arrival time within the 
frame t = mod T tends to be uniformly distributed over [0,T). Since 
the arrival process is memoryless, the unconditional density 


X. = lim P(X.- = j) 
n-H» '^n 

= P (arrival at steady state finds j packets in system) 
j = 0 , 1 , . . . , N 

can be obtained by averaging (28) and (29) over t, t uniform on [0, T). 
Thus the "actual" steady state packet system size density is 
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1 I 

X . J- ^ 

^ Xm I ° 1 - e 




K»0 


XTV^ „H ^ (XT)^ 


k“0 


i»jt+l 


N-1 

^■'-E 

j'O 


1 


j * 0, 1 N 1 (34) 


(35) 


where m Is given by (27), 

By averaging (30) and (33) over t we obtain the actual steady state 
message truncation and packet blocking probabilities respectfully; 


N 


vE^n.E 

i=o i=j+i 


(36) 


= <^)-'E^m e 

0=0 t=j+i 


(J^ - j)0o 


(37) 


The throughput s for the system is then 


s = (1 - Pb)p 


(38) 


The throughput s represents the average number of packet departures (or 
unblocked arrivals) per frame. 

At equilibrium, the average unblocked packet arrival rate, (1 - Fg) xF, 
equals the average packet departure rate, 1/m. Thus simpler expressions 
for the packet blocking probability and throughput are 
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( 39 ) 




)jSm 



(40) 


Before concluding this section, let us consider the "actual" packet system 
size when the messages consist of single packets. Evaluating (34) for this 
case we find: if g^ * 1, then 


X . « -4L j 

^ Am 


j * 0, 1, N - 1 


and 


This result could have been deduced from the following theorem [10]: for 
any stochastic system size process which changes only in unit steps, if 
either of 


ot- - lim P (unblocked arrival at t finds j in system) 
t*Ko 

ir.‘ - lim P (departure at t leaves j in system) 

^ t-Wo 

exists, then so does the other and they are equal. Note that the message 
system size process satisfies the conditions of the theorem (recall 
Figure 2b) even in the case of multi packet message arrivals. 

In the next section further utilization is made of equations (28) and (29). 
Both the packet and message steady state delay distributions are obtained 
by a straightforward application of these results, 
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IV. DELAY ANALYSIS 


The packet delay Is defined to be the total time the packet spends In the 
system* and similarly* the message delay is the time from its arrival until 
its last packet has completed transmission. This section presents results 
for the steady state packet and message delay distributions. Both the 
“virtual" and "actual" delay distributions are given. The term "virtual" 
signifies that the probabilities are conditioned on the arrival occurring 
at a time t s t mod T, while the "actual" distribution is unconditional. 

The results of this section require that an additional assumption be made 
on the STDM model defined in section II; the messages are transmitted in 
the same order that they are received in. That is, we assume a First^Come* 
First-Serve (FCFS) queueing discipline is used. 

Packet Delay; 

We will first find the steady state density of the virtual packet delay. 
This density will be representative of an "average" packet in exactly the 
same sense as discussed for the packet blocking probability. Thus, the 
order of service of packets in the same message can be assumed to be 
random (uniform). 

Let D(t) be the virtual packet delay for a packet which is contained in a 
message that arrives at a time A h t mod T, 0 < t < T. We set D(x) - +<« 
if the packet is blocked; the probability of this event is Pg(r) 

(eqn. (33)). 


The distribution of D(t) Is discrete with atoms at the points D(t) ® nT - t, 
n * 1, 2, N + I, and +«, We look for 

p(d(t) = nT - x|0(t) < -) = 


where 


P(D(t) <»)=!- Pg(t) 


(42) 


First consider the case 0 < t < T - 6T; 


For 1 < n < N 


^ P(X^ = j, L > n-j, (n-j)^^ served) 

P(D(t) = nT - t|D(t) < «) = r - '~ P ~ ( T) 

^ (43) 

where 

L = size of the message that the randomly chosen packet arrived in 


The density of L is given by (31). The probability that the packet will be 
served at a particular point in order ((n-j)^^) within a group of L is 1/L. 
Since the indicated events are Independent, 


n-1 


P^D(t) = nT - t(D(t) < «) = [1 - PgCx)] ’Ev" E 


7 > ^ 


(44) 


t=n-j 


0<t<T" 6T , l<n<N 

X, 
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Now consider the case T - 6T £t < T. In this case the delay now depends 
on whether the packets already in the system all arrived since the beginning 
of the current slot. If at least one arrived before the slot began, then 
a departure will occur at the end of the slot; otherwise a departure will 
not occur until a frame later. Thus we have for T - 5T < t < T 


P^D(t) = nT - t|D(t) < 


n-2 


[1 - ^ ^T-6T ' 2. n-j-1, (n-i-1)*^ served^ 

j=0 


n-1 


^T-6J > L > n-j, (n-j)^^ served)^ (45) 

j-1 


for 2 £ n £ N , while (ri = N + 1) 


p(d(t) = (N + 1)T - x|D(r) < 


N-1 


= [1 - Pb(t)]‘^ ^ P(x^ = j. = 0, L > N-j , served) 

(46) 


Since the arrival process is Poisson, it possesses the property of stationary 
independent increments, so that 

P(X^ = j, Xj_gj = 0) = Xq(T - 6T) p(j packet arrivals in [0, t - (T - 6T))) 
ji=0, 1,...,N , T-6 T£t<T (47) 
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where from (28) 


Xo(T - 5T) = 


V 


■XTO-6) 




(48) 


Separating independent events we may proceed as in the previous case, 
yielding: 

For T - 6T < t < T 


P^D(t) ” nT - t|D(t) < 


! n-2 «■ ^ 

j=0 A=n-j-l 


n-1 


^T-6T XI. 


0=1 


A=n-j 


(49) 


for 2 < n < N , while 


P^D(t) = (N + 1)T - t1D(t) < 


N-1 


= [1 - Pb(t)-' j]p(X^ = J. Xt.jt = 0) Z 


j=0 


t=N-j 


(50) 


where 


P(X, = J. X,.J, = 0) = x„(T - 6T) ^ e-xC-(T-^T)] (X[T . . , jT.z. . gm i.,*i 

i=0 


( 51 ) 
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and Xq(T - 6T) is given by (48). 

Together, (44), (49) and (50) give the steady state virtual packet delay 
density. 

Next we investigate the actual packet delay, which will be denoted simply 
D. As before we set D = +® if the packet is blocked. We wish to find 
the distribution 

P(D £ tglD < «) 

for all tQ. For each fixed tg, define (uniquely) ng and tq so that 

tg = IqT - Tg , n^ integer , 0 £ Tq < T (52) 

Let A be the time of arrival of the packet, and define t so that A = t mod T, 
0 £T < T. Conditioning on x we find that 

for T £ Tg 

P(D £ D < ®) = P^D(t) £ tg|D(x) < 

= p|d(t) £ ngT - t|D(t) < “I (53) 

while for t < xg 

P(D £ tg|x, D < “) = p|d(x) £ (ng - 1 )T - x1d(x) < (54) 
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Since T is uniform on [0, T), unconditioning on t gives 


P(D < tplD < ~) 


•'Tn 


t) < nJ - t|0(t) < “)dt 


if" P(«x) 1 

Jn 


(nQ - 1)T - t1d(t) < »)dT (55) 


or changing variables 


P(D <, tn|D < oo) = 


•^(nfl-DT 




(D(n(,T - t) iT - t|D(tl()T - t) < (56) 


This is the desired distribution function, and can be shown to be continuous 
even at the points t„ = n^T, n^ integer. It is differentiable except when 


tg = OqT; using Leibniz' rule we obtain 


P(D < tg|0 < ») = } P(0(To) = tglDUo) < -) 


(57) 


dt, 


for tg/T t integer , 


where t« is given by (52). 


Summarizing, the actual packet delay, conditioned on the packet not being 
blocked, has a steady state distribution given by (56). This distribution 
is atomless, and has a density given simply by (57). 


Message Delay : 

The message delay analysis is similar to the packet delay analysis, and is 
actually simpler since the order in which the packets are served is not a 
concern. Therefore only the definitions and final results are presented 
here. 

Let B(t) be the virtual message delay for a message Which arrives at time 
A = T mod T, 0 T < T. Set B(t) = +» if the message is completely blocked, 
thus, 


P^5(t) = +oo^ = Xhj(t) 


(58) 


The steady state conditional density of D(t) is given by (59)-(63): 


If 0 < T < T - 6T, 


n-1 


P^D(t) = nT - t|D(t) < = [1 - Xj^jIt)] ^ Xj(r) g^_j 

j=0 


for 1 < n < N - 1 , while 


(59) 


N-1 °° 

p(D(t) = NT - t|D{t) <») = [!- x^(t)]"'’^ Xj(t) ^ g^ 

j=0 Ji=N-j 


(60) 
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If T - 6T < T < T, 


P^D(t) = nT - t|D(t) < 


n-2 


= Cl - XN(t)r^|;^P{X^ -0, = 

j=0 


n-1 


■^Xj ^-6T “ %-j 

• •• 

J=1 


2 1 n < N - 1 , while 


p(D(t) = NT - t|D(t) < ») 


N-2 


= [1 ^ XN(t)]-^|^ P(X^ = 0, X^.,^ = 0) 
j=0 


(61) 


N-1 


*X1 - "‘’‘t = “• *T-6T = «)] E J 

j=l )l=N-j j 


and 


p(5(t) = (N + 1)T - t|D(t) < co) 

N-1 

- [1 - X|^(t)] ^T-6T " 9c 


j=0 


where P(X^ = j, X^_^^ = 0) is given by (51). 


i^=N-j 


(62) 


(63) 
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Let 5 be the actual message delay. The steady state distribution of 5 , 
conditioned on 5 < «, is atomless and has the density 

^ P(D < t^lD < «) = I p(d(tq) = tQ|D(TQ) < -) (64) 

for tg/T f integer, where Tq is given by (52). 

In this section we determined the virtual packet and message delay distribU" 
tions, expressed in terms of the virtual packet system size probabilities 
(28). The actual delay distributions were then given in terms of the 
virtual delay densities. 

If we were only interested in the mean actual packet delay E(D|D < «), it 
could be easily computed using Little's Result [11]: 

sE(DjD < ») = T E(X) (65) 

where s is the throughput (38) and E(X) is the mean of the actual packet 
system size density (34), (35). This simple result does not depend on any 
assumptions on the capacity N. Unfortunately, if N < », any higher moments 
of the actual packet delay must be computed by a numerical integration 
involving the density (57). In addition, a finite capacity constraint 
imposes a similar numerical burden in obtaining any of the actual message 
delay moments. These numerical difficulties are alleviated in case N = ®. 
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The next section concentrates on the Infinite capacity system. Transform 
expressions are found for the steady state system size and delay distribu- 
tions* These transforms are useful for determining moments; in particular, 
moments of the virtual delay distributions can be expressed in terms of those 
of the virtual packet system size distribution, 


V. TRANSFORM EXPRESSIONS FOR THE CASE N = » 


Further analysis is possible when the STDM system has an inffnfte buffer 
capacity, In this case simplifications occur in the transfomi representa- 
tions of the various probability distributions of interest. The transforms 
can be used to obtain moments; some results are given in this section. 

First consider the stationary distribution of the MC {Y^}. Define for 
W 5.1 

00 

g‘(z)=5Zv* 

to be the Z- transform of the number of packets per message. Also define for 
|Z1 < 1 



(67) 


and for i > 0 



( 68 ) 


The latter equality holds because is a function only of j-i (see (2)) for 

* J 

i >0. Using a straightforward generalization of Welch's result [12] to the 
case of compound Poisson arrivals, we find that the stationary distribution 
irr.) exists uniquely if p < 1. In this case the Z-transform of {tt.} is 

J w 

given by (|Zl i 1) 
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'(Z) 


j»0 


* ITn y 


y^(z) - zy^(z) 

'fb(z) - 2 


( 69 ) 


where 


_-aT(V-6) 


pe 


(70) 


Using (2) and (3) we find 


•k 

Ye(Z) 




1 * 


•XT(1-5)G*{Z) e-AT(l-G*(Z))]} 

^ 


and 


y;(Z) = 


(71) 


(72) 


so that 


_ (1 - o) e-^T«(1-G*(Z)) (■, _ ,-AT{l-G*(Z))l 
■ p-|I-if{i--G*iz)'; . 


(73) 


Now define the Z-transforrii of the steady state packet system size density 
at time t = t + ml (large integer m) to be ( |Z| < 1 , 0 £ t < T) 


'(Z. t) zJ 

j=0 


(74) 
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The mean steady state sojourn time Is from (27) and (70) 


m ^ T/p (75) 

Using this and (28) we find for p 1, |Z| < 1, 0 < t < T 

X*(Z, t) Ml - p)e"^<^^«T)0-G*(Z)) ^ py*(2),-Xt(l-6*(Z)) 


.a 


p) (1 - 2)e->^irHT)(UGHl)) 
■^T(T-6*(Z)) 


(76) 


Equation (76) provides a transform expression for the “virtual" packet 
system size probabilities at any time (during steady state conditions), 
and Is primary to all subsequent developments in this section. 

Note that X'*'(Z, T-6T) Is the Pollaczek-Khinchen transform equation for the 
system size at departure Instants of a bulk arrival M/D/1 queueing system. 

If we let 6 = 0, the Z-transforms of the packet system size just before and 
just after framing instants, X*(Z, T") and X*(Z, 0), agree with the results 
given by Konheim [3] and Hayes [4] respectively. 

Differentiating (76) with respect to Z and evaluating at Z = 1 shows that the 
mean packet system size is the linear function of t 

E(X^^mT^ = (1 - p)p6 + pT + Xt& (77) 


where 


1^1 

_ ~7 

where G and G are the first and second moments of the group size G. 


( 78 ) 
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The Z-transform of the packet system size density as viewed by an actual 
arrival at steady state Is defined as 



(79) 


and Is found using either (34) or by averaging (76) over t to be (p < 1 

|Z| < 1) 

[n - ^ PV.(Z)] 

(1 - o) n - fl - 


(80) 


The mean (actual) packet system size Is found to be 


E(X) = pS + p) 


(81) 


We now determine the transforms of the delay distributions (assuming FCFS), 
using the results of section IV. Define for 0 £ t < T, Re(s) £ 0 

D*(s, t) = Ete''“(^'] (82) 

to be the Laplace-Stieltjes Transfoim (LST) of the virtual packet delay 
steady state distribution. From (4S)-(50) and (70) we find that 

If 0 £ T < T - 5T 

D*(s, t) = 1 - G* . (e ^^) x*(e‘^^, t) (83) 

( 1 - e-®^) ^ 
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while if T - fit < T < T 


0*(s. t) “ e-s(T-T) |x*(e"®^, t) 


- (1 - p)0 - 

Rubin [7] pointed out that 


( 84 ) 


G*(Z) 


1 - G*(Z) 
(1 - Z) G 


(8B) 


is the Z-transfom) of the Backward Renewal Time (BRT) 6 of the renewal process 
with inter-renewal -times {G^}. For an arbitrarily chosen packet, this BRT 
simply represents the number of packets in the same message that are served 
in front of the chosen packet. 

Equations (83) - (85) allow moments of the virtual packet delay to be 
evaluated from moments of its components. For example, using 


E(G) = 


G^ - ^ 
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(86) 


and (77), we find 


E[D(t)] “ ^ “ T + XtIj - I* + T 


? 


Z(1 - p)G 


+ pT6 + 


0 if T < T - 6T 

(l-p)T if T > T - 6T 
(87) 


Now consider the LST of the actual packet delay, defined as 
D*(s) = ECe'®”] 


( 88 ) 


This can be determined from (83) and (84) by unconditioning on t: 


D.(s) = 0 - P)(l - e-^^) G«(e-^T) 
ST - XT [1 - G*(e"'^)] 


( 89 ) 


The mean actual packet delay is 


E(D) = 6T + T 


2(1 - p) G 


(90) 


Unfortunately? a simple relationship between D*(s) and X*(s) does not occur, 
except in the case of single packet messages. If g-| = 1, we obtain the usual 
M/6/1 result 


D*(x - XI) = X*(Z) 


(91) 


("generalized Little's Result"). 


Next define the LST of the virtual message delay as 


D*(s, t) = , 0 £ T < T 


(92) 


From (59) -(63) we find 


D*(s,t) = e^”^ G*(e"^^) X*(e“^^,x) 


1 


if T < T - 6T 


^-ST+AT[1-G*(e-=‘T)] ,f^,T-6T 

(93) 


Define the LST of the actual message delay to be 


D*(s) = ECe"®*^] 


(94) 
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Averaging (93) over t we find 


B.(s) = - P)(l - 

ST - [1 - G*(e-'T 


eSTQ-6) 

)3 


The mean actual message delay is 


(95) 


2 

E(D) = «T + TG - i + (96) 

^ 2(1 - p) G 

This concludes the analysis of the STDM/TDMA system. In the next section the 
various results presented will be illustrated by numerical examples. 
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VI. NUMERICAL EXAMPLES 


In this section several examples are presented in order to illustrate 
applications of the preceding theory. All cases refer to packet behavior 
at steady state. Examples of message behavior are coital ned in Van [8] 
and Lam [5]. 

Figure 4 gives an example of the steady state packet system size probabilities 
versus time for each possible state, computed from (28) and (29). A vertical 
slice taken at a particular time t = ml + x will yield the complete density 
{x.(t)} for that time. Also included on the graph is the virtual packet 
blocking probability Pg(T), represented as a dashed line. The nonstationary 
(periodic) nature of the system size process is clearly evident. The 
example used for Figure 4 is a source allotted 6 = .5 of each frame, with 
traffic Intensity p = .85 and a capacity constraint of N = 5 packets. The 
message length is fixed at 2 packets, so that gg = 1. Note that this causes 

x-j < )<2 < Xq. 

Figure 5 shows the mean system size versus traffic intensity for the simple 
case of single packet messages and an infinite capacity. The curves are 
parameterized by the slot-to-frame ratio 5, which determines the degree 
of "funneling" in the system. This figure indicates that as 5 is decreased, 
its incremental impact on the system behavior decreases rapidly. 

Figures 6, 7 and 8 were generated for a source allotted 5 = .01, with 
either gg = 1 or P 2 = 1 (i.e. packets arrive in either groups of 8 or in 
groups of 2 respectively). For each of these configurations the ;sy Stem 
capacity is either 10, 20, 40, or an infinite number of packets. 
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Figure 6 , like Figure 5, illustrates the mean ("actual") system size as a 
function of the traffic intensity p. We see that more congestion occurs 
for the more bursty input stream (gg = 1). The two cases cross over (same 
N) for large p values due to the relatively larger number of packets blocked 
in the case gg =1 (see next figure). 

The packet blocking probabilities Fg are given in Figure 7 as a function of 
p. Notice that the gg = 1 , N = 40 and gg = 1 , N = 10 cases coincide, con- 
firming the intuitive notion that a system with traffic four times as bursty 
requires four times the capacity to obtain the same Fg. 

Figure 8 shows the mean packet delays versus p. These curves were generated 
by dividing the mean system size (Figure 5) by the unblocked arrival rate 
(Little's Result); i.e., the mean is conditioned on D < ». 

In the next example we consider two input streams with the same mean number 
of packets/message F, but which differ in higher moments. Specifically, in 
one case g^ = 1 . while in the other we have gg = 1/2 and 84 = Og " 1/4- The 
mean (unblocked) packet delay versus the throughput s is given by Figure 9. 

A smaller delay is attained by the zero variance case, g 4 = 1 . Note that 
no crossover occurs for fixed N as it did in Figures 6 and 8 , because s is 
used instead of p for the abscissa. 

Figure 10 presents the actual packet and message delay densities for the 
same example as was used for Figure 4. Because of (57), this graph can 
also be used to determine the virtual delay densities. 
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V. SUMMARY 


An exact analysis of the STDM/TDMA system with Poisson message arrivals 
has been given. The results were obtained by rigorously defining the under- 
lying processes involved and then applying the theory of semi -regenerative 
processes. Steady state packet and message queueing behavior is predicted 
for possibly capacity-limited systems. 

The embedded chain is defined as the system size at departure instants. 

The transition probability matrix is given, from which the stationary 
distribution can be obtained using either Van's [8] method (N < «J or 
transform methods (N = <»>). The steady state system size at all times is 
then determined, arid is shown to be time-periodic. By averaging over the 
period, we obtain the system size density as viewed by a typical message 
or packet arrival. 

The system size densities are used to obtain the blocking and truncation 
probabilities (if appropriate). In addition, the delay distributions are 
expressed in terms of the packet system size density. 

Numerical examples illustrate the probabilistic behavior of the STDM/TDMA 
system. The periodic nature of the queueing process is exemplified, and 
its effect on the delay density. Other examples present the mean system 
sizes, mean delays and blocking probabilities for various system parameter 
values. 

The results of this paper will enable a performance prediction of the 
STDM/TDMA system. Considerable flexibility is allowed in modelling the 
packet arrival process. The design engineer can determine what buffer 
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size and allowable traffic load Is required to maintain given blocking 
and delay constraints. 
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APPENDIX A; MESSAGE SYSTEM SIZE ANALYSIS 


The STDM/TDMA model is essentially the same as for the packet system size 
analysis. The parameters T, 6, A, {g^} and p are defined as before. The 
remaining parameters are defined similarly except that they are measured in 
messages rather than packets. For example, 

N = capacity of the system in messages 
= system size in messages at time t 
= time of nth message departure 

etc. It is hoped that this duplication of notation does not cause 
confusion. 

Note that for N < ® the system inherently operates in a different manner 
than in the packet capacity-limited case. Messages that arrive to a full 
system are completely blocked, so that messages which are accepted into the 
system are transmitted in their entirety. If N = » the system operation is 
identical to that of the previous model, except now messages are counted 
instead of packets. 


The discrete process {Y^}, = Xp+ is a MC with transition probabilities 

(.see Yan [8]): 


i > 0, j < N - 1 : 






-UT 

(j-i+D! 


j - i + 1 > 0 


otherwi se 


(A.l) 
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j < N - 1: 


^Oj * 7-W 

i-1 ' " ® 

. [XTU-H«)]'**’)| 


(A.2) 


and for any i 


^iN-l 


= 1 


N-2 

E 

j=0 




(A.y) 


The MC {Yj^} Will possess a stationary distribution {ir^} if either N < « or 
p < 1, and will satisfy equations (6) and (7). 


The mean sojourn times are 



(A.4) 


and for j > 0 
m. = G T 

J 

so that 



(A.5) 


(A.6) 


Since is semi -regenerative, (21) and (22) will remain true. We now 
have for 0 < T < T: 
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(A.7) 


while for m > 0 (m Integer) 

Pq(X^ * j. Ri > t + ml) » P(G^ > m) 


+ ^ (jJ d -: t] )'^.. 


(A. 8) 


and for k > 0, m > 0 


/ -l(T+mT) p(G, > m) J - k 10 

^k^\+mT " ^ 

10 otherwise 

(A. 9) 

Combining (21), (22), (A.6)-(A.9) we find that the steady state message size 
density at time t = t + ml, 0 < t < T is 




g-AT(l-6) » g-x[(ji-l)T+6T+T] (x[(a - 1)T + 6T -f t])^ 

1 - ^ 

A=1 


i *kZ "‘“i * 


(A. 10) 


k=0 m=0 


0 < j < N, and 


N-1 

j=0 


(A.ll) 
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Of course, there is no message truncation, and X|^(t) represents the virtual 
message blocking probability. 

Averaging (A. 10) over t (uniform) we find after simplification 



But it can be easily shown that 
.i N-1 

° "i 

k=0 i=j 


(A.12) 


(A. 13) 


so that the steady state message system size density as viewed by arrivals 
is simply 


X = 

xi 


0 < j < N 


(A.14) 


with the actual message blocking probability 



j=0 


xm 


The system throughput is 

s = (,1 - 3L)p ~ ^ 

xm 


(A. 15) 


(A. 16) 
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The Case N « 

If the system has an Infinite buffer capacity* Lam L6] showed that (|Zl £l) 



s*(> - xz) - z s*(x - XI) 
Sj(j - M) - I 


(A.17) 


Where 


n rp) ( 1 ■ - e'^^) 
-AT(l-d) 


XT e 


(A. 18) 


and where 


sj(u) = G*(e-“^) 

S*(u) = r, . e-(>-“)(T-<ST) . ,-uT,1 

X-u 1 - J 


(A. 19) 
(A. 20) 


are the Laplace-Stieltjes transforms of the message "service times", 
conditioned on whether the message arrives to a busy or an empty system. 
Combining (A.17)-(A.20) yields 


Using (A. 10) we find (lZ| £ 1) 

X*(Z, t) = (1 - 

. ,t--XT(1-Z) 1 - 

+ ^Te Y (Z) -Y-^^t(i:z) ' 


(A. 21) 


(A. 22) 
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The mean steady state message system size at time t - t + mT, m Integer, 

0 T < T Is 


W> ■ »(t - T + ST) + ptV + T - p + 1 T (1 - s)] + - UT)^ I? ■ g) 


where 


(A. 23) 


T - p - + »TS + 

By summing (A.14) over all j we find that Am « 1 and that 
Xj » TTj for all j , 


(A. 24) 


or 


t{Z) = Y*(Z) 


(A. 25) 
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APPENDIX B: DISTRIBUTION OF ARRIVAL TIME MOD T 


It Is Intuitively clear that the limiting distribution of the arrival time 
within a frame Is uniform when the arrival process Is Poisson. A formal 
proof of this fact Is given here. 

As before we let 

A^ = time of the n^^ arrival 

and assume that arrivals form a Poisson process of rate x. Let 

I— " A_ ■■ A_ I 
n n n-1 

be the interarrival times, and let 

'n ' - TLVTJ 

be the time from the last frame Instant to A^^, where T Is the frame duration 
and notationally 

LxJ = greatest integer less than or equal to x 

It Is clear that the sequence of random variables defined on the state 
space [o,T), forms a Markov Chain (MC). We compute the transition probabilities 

P(Vl ^ ^1 I "^n ' ^0^ 

by considering two separate cases. 
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P(v, < t, 1 - t^) - P(0 i - t^) 


".4^ " ^0 -< ^n+l "’T “ ^0 " 

m®l 


i’) ^}<^o 


= 1 - ■ ^0^ 


j“>(to ■■ t] ) 


1 - e 


2 ^1 i = ^0^ ^Yj " ^0 ^n+i ^ “ ^0 

m=l 


g"^(T ~ tp) j-1 _ Q~^t~|' 


1 - e 


Given the distribution of the distribution or can be found via 


*'^Vl - *1 


> = /pl Vl 2 ‘l I 'r = to> £ *o> 


Suppose is uniformly distributed on [0,T). Then for t-j e [0,T) 

’’<vi£*i>= / ' n - ^ 

Jq 

, e-MT - to) [, . j-xt,j dt„ 

L ^ - e-^T T 


T 

so that is also uniform on [0,T). Thus the uniform distribution is 
stationary for the MC. By Proposition 3.2.10 of [18] we know the stationary 
distribution is unique, and by Orey's Theorem (Theorem 6.2.8 of [18]) we know 
the distribution of the MC converges to it. 



Figure 1. STDM Structure 
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Figure 2b. A Message System Size Sample Function, Finite Message Capacity 



Figure 3. The Semi-Markov Process 
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Figure 4. Packet System Size Density vs. Time 



Figure 5. Mean System Size vs. Traffic Intensity for 
Various Normalized Slot Widths 
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Figure 6. Mean System Size vs. Traffic Intensity for 
Various System Size Capacities 
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Figure 7. Packet Blocking Prohab.il ity vs, Traffic Intensity 
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Figure 9. Mean Packet Delay (in units of frames) vs, Throughput 
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